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, Quantum Monte Carlo simulations were carried out for the system of a magnetic trimer on a metal 

surface. The magnetic trimer is arranged in two geometric configurations, viz., isosceles and equilat- 
eral triangles. The calculated spectral density and magnetic susceptibility show the existence of two 
regimes: Yosida-Kondo dominant regime and magnetic frustration dominant regime. Furthermore, a 
critical crossover between these two regimes can be induced by changing the configuration of the 
magnetic trimers from isosceles to equilateral triangle. 
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S 

Q i Studies on systems of magnetic monomers and dimers on a metal surface have allowed us to gain 

deeper understanding of, and also clarify, the nature of typical phenomena associated with the classical 



'sj" \ Kondo systems, e.g., the real-space image of Yosida-Kondo resonance and the competition between 

0> ■ the spin-flip interaction and the indirect Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction. ° Fur- 

\C) , thermore, in the system of a magnetic trimer on a metal surface, we found that changing the geometric 



configuration of the trimer, from isosceles to equilateral triangle, changes the corresponding spectral 
density near the Fermi level. In the isosceles configuration, the spectral density exhibits a sharp peak 
near the Fermi level, which we attribute to the Yosida-Kondo resonance, with a corresponding rela- 
tively high Kondo temperature. While in the equilateral configuration, no peak is observed near the 



^ ■ Fermi level.^ These results are in agreement with scanning tunneling spectroscopy (STS) experiments 

on Cr trimer on Au(lll).^ These results also indicate the existence of two separate regimes. One is 
a Yosida-Kondo dominant regime for an isosceles trimer where the ground-state of the system is sin- 
glet. The other regime is a magnetic frustration dominant regime for an equilateral trimer, in which 
the ground-state of the system is degenerate due to the competing antiferromagnetic (AF) interactions 
between adatom pairs. 

Can we observe the crossover between the Yosida-Kondo dominant regime and the magnetic frus- 
tration dominant regime? This is the main concern in this paper. To answer this question, we again 
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consider the system of magnetic trimer on a metal surface, as shown in Fig. 1. By calculating the 
spectral density, as well as the magnetic susceptibility at adatom 3 for a wide range of temperature, 
we clarify the existence of two separate regimes and suggest the nature of the crossover between these 
two regimes. 



Fig. 1. Magnetic trimer on a metal surface. Separation between adatoms 1 and 2, Ru - a. a is the surface 
lattice constant. Adatom 3 is moved along the direction of the arrow, maintaining the relation R13 - /?23- 



A system of three magnetic adatoms on a metal surface is modeled as follows, 

i i cr,(ij) 

+ ^ fk^ko- + Y^iVkicl^dia- + h.c). (1) 

ko" k(T,i 

Adatom indices /, j = 1, 2, 3. Spin index cr =1, J,, e; is the energy level associated with the electrons at 
adatom site /. Ui gives the Coulomb repulsion at individual sites, tjj is the direct hopping between the 
adatoms (/ j), and gives rise to the direct AF interactions, k is the metal surface conduction electron 
wave- vector. Vkr = e '^^' Vq gives the hybridization between the states localized at the adatoms and the 
conduction electron medium. R, is the adatom position. We choose a free-electron dispersion relation 
for the 2D conduction electron medium, i.e., 

D k , 

ek= 2-[(^)'-l]. (2) 

D is the conduction electron bandwidth and is the Fermi wave-number. Since the wave-number is 
within the range < ^ < V2^f> then - ^ < eu < ^. 

The Green's functions are defined so as to satisfy the equation 

^(/6J„ - H)i,yGyi,{iaj„) ^ 6^1,, (3) 

V 

with adatom indices /, j and wave- vector of conduction electrons k given by ji, v, and k. Un - {2n + 
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l)7r/y6. n is an integral number. /? - l/k^T (eV)"^ is the Boltzmann constant. From eqn. (3), we 
derive the non-interacting Green's function, which corresponds to Uj = case, and get 



G°.(/w„) = 



iaj„ - 61 - f -hi - F12 -tis - -F13 

-til - F2I i(^n - ^2 - Fo -t23 - F23 
-t31 - Fsi -t32 - f 32 i(^n - h - Fq 

with e,- = 6; + Ui/2. Fq and Fij (i j) are as given in our previous study. ^ 

AD xdx 



(4) 



^ lix)n - ^ Jo 



ioJn - ^{X^ - I) 

Vle-'^^" AD xJoixkpRiAdx 



^ _y Vp-e— V _ AD r 



/W„ - -yCx^ - 1) 



(5) 



(6) 



R,-y = R,- - Kj, A = V^kp/ID, and /o(-^^f^12) is the order Bessel function of the first kind. 
Fo{ia)„/A) accounts for the finite peak width at e,-, e,- + Uj, and also the Yosida-Kondo resonance at 
the Fermi level, at low temperatures. Fij{ia)n/ A,kpRij) gives rise to the RKKY interaction between 
adatoms. 

In this study, we used quantum Monte Carlo (QMC) simulations to get the interacting Green's 
function corresponding to the Uj + case and the magnetic susceptibility at each adatom.'""'^ We 
used the maximum entropy method to get the spectral densities from the Green's functions. 

From hereon, we consider the specific case of a Cr trimer on Au( 111), in the configuration depicted 
in Fig. 1. The separation between adatoms 1 and 2 is kept constant, while adatom 3 is shifted. We 
assume that the direct hopping between adatoms ?,y is inversely proportional to adatom separation Riy 
Since adatoms 1 and 2 are fixed, we can set t\2 equal to a constant. By further setting = ?23 - 
fi2 x7?i2//?i3, the direct AF interaction rising from can be tuned by varying R^j. The corresponding 
spectral densities at adatom 3 are shown in Fig. 2. 

When 7?i3 = R23 » R12, adatom 3 on Au(lll) forms a single impurity Kondo system with 
a corresponding extremely low Kondo temperature. At high temperatures, the magnetic moment of 
adatom 3 is not fully compensated by the conduction electrons, therefore no Yosida-Kondo resonance 
would be observed.''^ The dimer, which consists of adatoms 1 and 2, is in the singlet state {H} due to 
the strong AF interaction. 

When adatom 3 is moved closer to the dimer, we found that the resonance at the Fermi level 
is most pronounced for the isosceles trimer with /?i3 = 5.0A, corresponding to ?i3 = 0.06eV (the 
red/solid curve in Fig. 2). We expect that the dominance of the Kondo effect at adatom 3 is enhanced 
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Fig. 2. (Color online) Spectral densities at adatom 3 in the system of Cr trimer on Au(lll) at 54K. The 
Fermi level is set as the energy reference. Parameters corresponding to Au(l 1 1) are set as follows: lattice 
constant a = 2.9 A, Fermi wave-number kp = 0.2A and the half-bandwidth D/2 - O.SeV.'"* For Cr 
adatoms, the energy level of d electrons ed = -O.lSeV, Coulomb repulsion U = 0.3eV, and tn = O.lOBeV. 
A = O.OlSeV results from the coupling between adatoms and surface and gives the broadening of d levels.^ 
As kpRi2 = 0.58 and ^f^13 < 1 .6 then the RKKY interaction between adatoms rising from F,y in eqn. (4) is 
ferromagnetic.^ In the system specified above, the magnitude of the AF interaction rising from is much 
larger than that of the RKKY interaction, then the calculated densities are shown as f 13 dependence. 

by the combination of AF interactions between adatoms, since adatom 3 is dependent not only on the 
magnetic interactions with other adatoms, but also on the interaction between adatoms 1 and 2 J The 
ground state of this isosceles trimer system is {f J, •!, where {•) stands for the Yosida-Kondo singlet. 

When /?i3 = R\2 corresponding to t\j, = O.lOSeV, the high peak structure at the Fermi level 
disappears (the blue/dot-dashed curve in Fig. 2). The competition of the equally strong AF interactions 
between the adatom pairs gives rise to the dominance of magnetic frustration. In terms of the spin 
model, H = Z<,j> JS jS j with 7 > and S = +1/2, the ground state of the equilateral trimer can 
either be (TJ-t) or {^ii), associated with the energy E = -J/4 lower than the Fermi level, while the 
excited state {tttl associates with the energy E = 37/4. Therefore, we observe the Yosida-Kondo 
resonance suppressed due to the dominance of magnetic frustration in the equilateral trimer. 

In Fig. 2, as ?i3 increases from 0.06eV, the spectral density at the Fermi level drastically reduces. 
We observe the resonance feature shifting away from the Fermi level and merging with the d reso- 
nance. The small peak rising near the Fermi level is the Yosida-Kondo resonance. This behavior is 
different from that observed in the dimer system, where the Yosida-Kondo resonance splits into two 
nearly symmetric (with respect to the Fermi level) peaks. '^'^ Thus, we expect a difference in the nature 
of the crossover behavior of a trimer system and a dimer system. In the dimer system, the crossover is 




4/8 



J. Phys. Soc. Jpn. 



Letter 




low T 



'1.2 



- 1 



■ 0.8 



■ 0.6 



0.4 



0.2 



2x10" 



5x10"^ 1x10° 2x10° 



5x10" 1x10 



Fig. 3. (Color online) Contour plot of the spectral densities at the Fermi level at adatom 3 versus temperature 
and magnetic interaction, the contour spacing is 0.14. At each blue (white) point, the spectral density as a 
function of energy is presented in Fig. 2. 



from the Yosida-Kondo singlet state at each adatom to the AF singlet state {H). In the trimer system, 
the crossover at adatom 3 is from the Yosida-Kondo singlet state to the degenerate states as a result of 
magnetic frustration, in which spin direction can either be {tl or {J,}. 

In order to track the two separate regimes and observe the crossover between them, we study the 
temperature dependence of the spectral density at the Fermi level. Using the standard definition of 
AF interaction rising from the direct hopping t and scaling with A,'^ we have /12 = 4t^2/^^^ ^^'^ 
Ju - 4fj3/?7A. From this point on, we use //12 as a variable, instead of fi3, in order to present the 
dependence of properties at adatom 3 not only on /13, but also on the interaction between adatom 1 
and 2. The equilateral trimer corresponds to = Jn ~ 10, while J^^/Jn ^ 1 corresponds to 

the case when adatom 3 is sufficiently far from adatoms 1 and 2, and forms a single impurity Kondo 
system with k^T^"^'^ /A = n x 10"^. Within the intermediate range T^j/Jn, we show the spectral 
density at the Fermi level in a broad range of k^T/A in Fig. 3. The spectral densities at the Fermi level 
for various ?i3 shown in Fig. 2 correspond to the blue (white) dots in Fig. 3 at k^T/A w 0.3. 

When //12 > 2, the corresponding density at the Fermi level is almost zero, not only at a finite 
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Fig. 4. Temperature dependence of the magnetic susceptibilities at adatom 3 (solid lines) at magnitudes of 
interactions corresponding to the arrows indicated in Fig. 3. Dashed lines depict the universal shape of 
Kondo susceptibility with the effective Kondo temperatures k^T^A = 0.1 and 0.2.'^ '^ The error at each 
calculated point is smaller than or equal to the radius of the solid circle in the figure. 



temperature k^T jlS. w 0.3, but also in a broad range of temperatures. In the range 2 > fl-^/Jn > 
1, the lower the temperature, the more significant the increase in density with decreasing J^^/Jn- 
When J^^/Ju < 1, the density decreases more gradually at lower temperatures. Thus, we observe two 
separate regimes, viz., the magnetic frustration dominant regime in J\^IJn ^ 2 where the density at 
the Fermi level is nearly equal to zero; and the Yosida-Kondo dominant regime in the range J^^^IJn < 
1 where the density at the Fermi level increases as the temperature decreases. The crossover point 
between two regimes should appear in the range 2 > J^t^/Ju > 1 where the variation of density versus 
J^^/Ju is drastic at low temperatures. 

To further clarify the nature of the crossover between the Yosida-Kondo dominant regime and 
the magnetic frustration regime, we study the temperature dependence of magnetic susceptibility at 
adatom 3. The magnetic susceptibility at adatom 3 is defined as ;\^3 = dT{cr^{T)o-j(0)}}^ (Ti = 
- n^.. Fig. 4 shows k^Tx^ versus k^T/A for various Jf^/Jn, corresponding to the aixows in Fig. 3, 
with the corresponding order from right to left. These plots exhibit three different behaviors, indicating 
the existence of two separate regimes and the crossover between them. 

JlT,IJn ^ 2 corresponds to the magnetic frustration dominant regime. The corresponding curves 
of k^Txs with decreasing k^T/A do not follow the universal curve of Kondo susceptibility, but de- 
crease rapidly and smoothly. 

J^^/Jn ~ 1 belongs to the Yosida-Kondo dominant regime. The corresponding curve of k^Txi 
versus k^T/A follows the universal curve of Kondo susceptibility for a single impurity case with the 
effective Kondo temperature k^T^/A - 0.1.'^'^^ 
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k^Txi versus k^T/A, corresponding to 2 > j\-^IJu > 1, also follows the universal curve with 
k^T^/A = 0.2 at high temperatures. However, as T < T^, k^Txs versus k^T/A does not continue 
following the universal cui"ve, but instant drops rapidly close to zero. 

In summary, as we move adatom 3 towards the AF quenched singlet of adatoms 1 and 2 so as 
J^^/Ju < 1, adatom 3 coupling to the metal surface effectively behaves as a single impurity Kondo 
system but with a higher Kondo temperature than that corresponding to //12 1. In the particular 
case of J^^/Ji2 ~ 1, k^T^/A = 0.1 » k^T^"^'" /A = n x 10~^. And as expected, the crossover is 
found to occur in the range 1 < J\^IJn < 2 via the study on the temperature dependence of k^Tx?,- 
At a finite value of J^^/Ju, we observe the discontinuous plot of k^Txs versus k^T/A. Therefore, we 
suggest that the critical crossover would be observed in a system of a trimer on a metal surface. As 
10 > //12 > 2, no Kondo temperature is defined due to the dominance of the magnetic frustration. 

In conclusion, we found the existence of two separate regimes in a system of a magnetic trimer 
on a metal surface: the Yosida-Kondo dominant regime and the magnetic frustration dominant regime. 
Moreover, we suggest the critical crossover between these two above regimes. 
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